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Abstract 

Spectral Spatial Fluctuations (SSF) of the Cosmic Microwave Background Radiation 
(CMBR) temperature are considered as a result of an interaction of primordial 
atoms and molecules with CMBR in proto-objects moving with peculiar velocities 
relative to the CMBR. Expected optimistic values of AT/T are 2 x 1CT 5 -2 x 1(T 6 for 
SSF caused by HeH + at z =20-30 which are possible redshifts of early reionization 
scenario. The bandwidth of the lines is 0.1-2% depending on the scale of proto- 
objects and redshifts. For the SSF search CMBR maps in different spectral channels 
are to be observed and then processed by the Difference method. Simulation of the 
experiment is made for MSRT (Tuorla Observatory, Finland) equipped with a 7 x 4 
beam cryo-microbolometer array with a chopping flat and frequency multiplexer 
providing up to 7 spectral channels in each beam (88-100 GHz). Expected AT/T 
limit in the experiment is 2x10 ~ 5 with 6'- 7' angular and 2% frequency resolution. 
Simulation shows that SSF may be recognized in the angular power spectrum when 
S/N in single frequency CMBR maps is as small as 1.17 or even something less for 
white noise. Such an experiment gives us a possibility to set upper limit of SSF in 
MM band and prepare future SSF observations. 
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1 Introduction 



Observational effects caused by primordial molecules seem to be most promis- 
ing in investigating Dark Ages epoch of the Universe. The basic properties of 
the molecules are discrete narrow energy levels and high efficiency of their 
interaction with CMBR. This leads to forming SSF if proto-objects, con- 
taining these molecules, move with peculiar velocities V p relative to CMBR 
(Dubrovich,1977,1982, Puy,1993). We may consider SSF as manifestation of 
the proto-objects at high redshifts 300 > z > 10 as scattering of CMB pho- 
tons by primordial molecules must leave imprints on the CMBR temperature 
distribution. 

The most abundant chemical elements predicted by the pure Big Bang model 
are H, He, D, 3 He, Li and their ions. They give some molecules in the primor- 
dial matter at z=100-200 such as H 2 ,H^,HD, HD+, HeH+, LiH, LiH+, H 2 D+. 
There are other molecules which should be considered if non-standard nuclear 
synthesis at the early times (z=10 9 ) or star formation at z = 100-200 took 
place. 

For the high efficiency of interaction between molecules and CMBR photons 
two values are important - the cross-section for scattering and the concen- 
tration (or relative abundance) of this molecule. The first parameter depends 
on the specific quantum structure of the molecule - its symmetry and dipole 
moment. The second one depends on the abundance of the chemical elements 
of which it is composed and on the rate of the appropriate chemical reac- 
tions. According to these constraints, we should take into consideration only 
those molecules which have a large enough dipole moment and relatively high 
abundance. These are HD + , HeH + , LiH, H 2 D + . We consider them as the 
basic molecules interacting with CMBR at high z (Dubrovich, 1977, 1982, Puy 
et al.,1993,Stancil et al.,1996). 



2 Theoretical predictions of SSF 

Let us remind physical conditions in Dark Ages epoch when the Universe was 
younger than 300 mil. years. The Universe was in (z + l) 3 / 2 times smaller and 
in (z + l) 3 times more dense than now, temperature of CMBR and molecular 
gas followed T z=0 (z + 1) law till z ^150. Then, temperature of matter should 
be lower than CMBR temperature due to different adiabatic indexes of ra- 
diation and matter (Varshalovich, Hhersonskii, Sunyaev, 1981) and it must 

Email addresses: dubr@MD1381.spb.edu (Viktor Dubrovich), 
bajkova@gao.spb.ru (Anisa Bajkova), vkh@brown.nord.nw.ru (V.B.Khaikin). 



2 



Table 1 



Temperatures of CMBR and molecular gas in Dark Ages epoch z=300-10 



z 


Temperature of 


Temperature of 


Commentary 




CMBR, 


molecular gas, 






K 


K 




1100 


3000 


3000 


Recombination 


300 


815 


815 


Beginning of Dark Ages 


150 


415 


415 


Middle of Dark Ages 


30 


85 


17 


Dark Ages + reionization 


20 


57 


8 


Dark Ages + reionization 


10 


30 


100-10000 


End of Dark Ages, first stars and quazars 



follow T z= o(z + l) 2 law till first star formation. Calculated temperatures of 
CMBR and molecular gas in Dark Ages epoch are given in Table 1, we used 
T 2=0 = 2.726 K (Mather et al., 1994). In epoch before secondary ionization 
molecules are formed in accordance with kinetic rate of specific chemical re- 
actions when matter and radiation are cooling. After reionization molecular 
formation follows a more complicated law due to growth of ionization degree 
of hydrogen and existence of heavy chemical elements produced by first stars 
(Basu, Hernandez-Monteagudo, Sunyaev, 2004). 

There are two mechanisms responsible for SSF formation: simple scattering 
(Dubrovich,1977) and luminescence (Dubrovich,1997). In the simple mecha- 
nism opacity in narrow molecular lines and peculiar motion of the proto-clouds 
results in CMBR disturbances. This can be more visible in the rest frame of 
the proto-object. The CMBR in this frame becomes non-isotropic and out of 
thermal equilibrium. On the side towards which the proto-cloud moves, the 
temperature of the CMBR will be higher than average and on the opposite 
side it will be lower. After reflection, the photons are distributed isotropically 
in this frame that leads to non-isotropic distribution in laboratory frame. This 
explains the principal role of opacity and peculiar velocity. The amplitude of 
the effect depends on spectral index of the reflecting radiation (Dubrovich, 
1977, Maoli et al., 1994). This effect corresponds to the elastic scattering be- 
tween the molecules and photons, i.e. the total number of the photons does 
not change. 

In fact all molecules have quite a complicated energy level structure. This 
allows for the possibility of a non-elastic process. It is the well known lumi- 
nescence process which plays an important role in the formation of radiation 
from a reflection nebula. Luminescence may enhance the intensity of rotation 
lines due to breaking of vibrational line quants (Dubrovich, 1997). 
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The magnitude of the CMBR temperature fluctuations due to the pure re- 
flection of photons by a moving object depends on peculiar velocity V p of a 
proto-object and optical depth r(Sunyaev and Zel'dovich, 1970, Dubrovich, 
1977, Maoli et al., 1994). The value of the distortions produced by molecules 
in CMBR is: 

AT/T = (V p /c) • t ■ K, (1) 



where c is the speed of light and r is optical depth of the proto-object ( r << 1, 
Vp/c < 10~ 3 ), K - gain coefficient due to luminescence effect, K=1-1000(K=1 
for simple scattering). For the linear stage of perturbation evolution V p could 
be estimated as (Zel'dovich, Novikov, 1975): 

V p = V p0 /(z+l) 1 / 2 (2) 



where V^ =600 km/s at large scales of CMB anisotropy (CMBA) and factor 
2-5 higher at small CMBA scales. 

Luminescence effect is tempting for SSF search with limited sensitivity but 
high K may be expected in narrow z intervals (Dubrovich, 1997) while the 
simple mechanism of SSF formation operates in the wide field of z. 

The optical depth r of the proto-object may be calculated as (Dubrovich, 
Lipovka, 1995): 

T = M 5/2 d 2 (l + zf' 2 n H a ■ [1 - exp~ h ^ kT }/3hH , (3) 

hv/kT = hc/[kT X (l + z)\ (4) 



H - Hubble constant, h-Planck constant, k-Bolcman constant, A - the rest 
wavelength, T - CMBR temperature, n# = n 0H ■ (3, where nn - density of H, 
n H - critical density of atomic hydrogen HI(10~ 5 ), /3 - a real density of H in 
proto-object which may be less or much more than unity (/3 > 4 • 10~ 2 ), d - 
dipole moment, a - abundance, a = 10~ 8 — 10~ 17 for different molecules and 
conditions in Early Universe. 

Let us consider LiH, H 2 D + , HeH + HD + primordial molecules. 

LiH is a very important molecule, because it consists of primordial Li. Its 
abundance is a good test for the epoch of nuclear synthesis in the Early Uni- 
verse. Its large dipole moment and relatively low frequency of the rotational 
and rovibrational transitions lead to a high value of K. But, unfortunately, 
some difficulties with the chemical processes of forming this molecule lead to 
a not too optimistic prediction for its abundance a = 1CT 11 — 1CT 17 . The rest 
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wavelengths of LiH are A r =676 //m(rotational transition), X v —7.1 /im (rovi- 
brational transition). 

HD + is also an important molecule due to the presence of primordial deu- 
terium D. The abundance of D is about 5 orders of magnitude larger than 
that of Li. But HD + has a dipole moment about 10 times less than LiH and 
a cross-section which is 100 times smaller. Another small factor is the abun- 
dance of H + at redshift z = 200, which might be about relative to that of 
neutral hydrogen. Due to the relatively high frequency of the rotational and 
rovibration transitions, the resulting value of K is not very large. Yet, if high 
sensitivity were reached, this molecule might be seen. The rest wavelengths of 
HD + are A r =227 /im, \ v =3.0 /jm. 

HeH + does not have any low abundance species. There are only two small 
factors which lead to a low abundance: a high rate coefficient for destruction 
(by electron recombination and collisions with the neutral atoms of hydro- 
gen) compared with the very small rate coefficient of formation, and small 
abundance of H + at high redshift. It might be the most likely molecule to be 
searched for due to more optimistic predicted abundance a = 10~ 8 — 10~ 15 . 
The rest wavelengths of HeH + are A r =149 /im, A^=3.3 jum. 

H2 D + is the simplest triatomic molecule with a high dipole moment. It con- 
tains primordial deuterium D. Due to the presence of very low frequency 
transitions in its spectrum the value of K can be very high. It is very impor- 
tant that the redshift of the H 2 D + recombination is relatively high. The rest 
wavelengths of H 2 D + are A r =1920 /im, \ v =4.5 /im. 

The expected abundances of these molecules in the Early Universe are dis- 
cussed by many authors (Lepp and Shull,1984, Puy et al, 1993, Palla et al, 
1995, Maoli et al, 1996, Stancil et al, 1996). 

One of the main observational effects caused by expansion of the Universe is 
the redshift of photons, expressed as \ Q b s = (z + 1) A em , where \ em is the 
emitted or rest wavelength of the photon. So, if an object contains a gas of any 
molecules and has a z redshift , it can only be seen at different but discrete 
frequencies z/j = u oi /(l + z), where v oi is a discrete set of the molecule's rest 
transition frequencies. And vice versa, if some object is seen at given frequency 
Vi this object can have the redshifts Zj = u oi /ui -1. 

The final spectrum depends on the angular size of the object 6. The angular 
size 6 and the frequency bandwidth Ais/is depend on linear size L (for a 
spherical shape object) and z of the object to be observed. For small z the 
angular size of the object (if the linear size is constant and equal to L) is 
decreasing with moving away, but it grows for z larger than some value. Exact 
relations between 9 and z in this model contain the Q m and Q\ parameters, 
which are the ratios of the matter n m and "vacuum" densities to the critical 
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density n c respectively (Q m + fl\ = 1). The fact that the cloud occupies the 
interval of redshifts z means that if it radiates or reflects the radiation locally in 
sufficiently narrow lines, all the radiation occupies the frequency interval Au, 
moreover Avjv— Az/(1 + z). Scattering by molecules in several lines leads to 
superposition of images of different proto-objects at different z. Taking into 
account Az dependence on L and O (Sahni, Starobinsky, 2000) we have: 

Avjv « 9<j>(z)[n m (l + zf + Ov] 1/2 ](l + z)~\ (5) 

2 

<P(z) = J + z) 3 + n A dz, (6) 



9 = (H L/c)(l + z)/(j>(z), (7) 
Avjv » (H L/c) yJn m (l + z) 3 + n A . (8) 



If we simultaneously observe CMBR in two receiver channels with frequency 
difference Au, the SSF with an angular size larger than some critical 6 C are 
seen at both receiver channels. If the size of the fluctuations is smaller, they are 
only seen in one of the receiver channels, another one will detect background 
noise only (Dubrovich,1982). In this case the correlation function of these two 
observations will be within noise to zero at small scales and different from 
zero at angular scales greater than 9 C . Then the difference signal of the two 
channels is not zero at the scales less than 9 C (L C ). Frequency separation Avjv 
of receiver channels determines a critical scale of proto-objects L c and scales 
L < L c remain in difference of two channels. Frequency resolution or the 
channel bandwidth Avjv determines the optimal scale L opt which has the 
best S/N at given z. The channel bandwidth should be a few times less than 
frequency separation and L opt < L c . 

The above consideration is correct for objects freely expanding over the Uni- 
verse in accordance with the Hubble law. But at the sufficiently late stages 
the deceleration of their expansion takes place and contraction begins due to 
self-gravitation. In this case the width of the line can be considerably smaller, 
and the amplitude larger (Zel'dovich, 1978). 

Fig.l. shows expected r(a, n#,z), AT/T(a,%,z), (fi(Q m ,Q L , z), 9(L,z), 
Av/v(L, z) for HeH + (A =1.49 10~ 2 cm, d=1.66 deb ) and K=l. According to 
Fig.l an optimistic level of AT /T = 2 x 10~ 5 — 2 x 10~ 6 may be achieved with 
a = 10" 9 — 10~ 10 . Expected observational frequencies f(z) are given in Fig.l. 
for rotational transitions(l=0) of LiH, H 2 D + , HeH + HD + . For (f>(fl m , £Il,z) 
calculation we used Q m and Ql which are discussed in several Cosmological 
models and recent CMBR experiments (Ruhl et al, 2004, Miller, 2004). 

Let us limit possible physical proto-object scales L in the Early Universe: 2 
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Fig. 1. ct>{n m ,n L ,z), r(a,n H ,z)6(L,z), Av/v(L,z) AT/T, f(z) for HeH+ 
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Mpc > L >0.002 Mpc. The most of Cosmological molecules considered may 
exist at z<200 and some at z<300. Then possible angular scales 6 for different 
Au/v>, L and z are calculated with (7), (8) and presented in Table 2. Table 2 
shows that modern angular scales 6 = 5' — W are still reasonable to search for 
large scale SSF at z=20-30 with Au/u=l-2%. Possible observational frequen- 
cies f(z) for HeH + (I=0) are given in Table 3. Most of primordial molecules can 
radiate at several rotational transitions I that leads to SSF formation at set of 
frequencies fi(z). LiH has radiating transitions at 1=0,1,2,3,4,5,6,7. The opti- 
cal depth r in molecular lines depends on I and in the case of thermodynamic 
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Table 2 



Angular and spatial scales for HeH + (I=0) 



Ai//i/,% 


Zmin 


Zmax 


L max ,Mpc 


-^min ,MpC 


@max ( ) 


@min( ) 


0.1 


5 


200 


0.5 


0.0025 


1.4 


0.14 


0.5 


10 


200 


1 


0.015 


4.3 


0.85 


1 


10 


200 


2 


0.025 


8.6 


1.4 


2 


15 


200 


2 


0.05 


11.5 


0.3 


4 


25 


200 


2 


0.1 


17.3 


5.7 


5 


30 


200 


2 


0.15 


20.1 


8.5 


6 


35 


200 


2 


0.17 


22.9 


9.6 


7 


40 


200 


2 


0.2 


25.8 


11.3 


10 


50 


200 


2 


0.25 


31.4 


14.2 



Table 3 

Observational frequencies f(z)in GHz for HeH + (I=0) 



z 5 10 15 


19 


20 22 


52 


60 


65 


f 335.57 183.30 125.84 


100.67 


95.88 87.54 


37.99 


33.00 


30.51 


Table 4 

Possible observational frequencies f(z) in 


GHz for LiH(L 


=6) 






z 10 15 20 


30 


32 34 


35 


80 


101 


f 287.08 197.37 150.38 


101.87 


98.68 90.23 


87.71 


38.9 


30.96 



equilibrium for optically thin layer is (Dubrovich,1982): 



r/,7+1 = T • (I + 1) • exp(-T I(I + 1)/T) (9) 

where optical depth t =t(I=0) is given by (3), T - temperature correspondent 
to energy of transition 1=0-1, T- temperature of radiation at given z. For LiH 
extremum of r lies at 1=6 (Ao = 0.95- 10~ 2 cm) and r(J = 6)/ro=2-4. Possible 
observational frequencies fi(z) for LiH(I=6) are given in Table 4. Expected 
AT/T level for SSF caused by LiH with existing estimates of LiH abundance 
is less than 10~ 5 that can not be achieved in this experiment but we can set 
an upper limit of AT/T for LiH at z=30 as well. 

If we realize a channel bandwidth of Av/v—1% in 3 mm band preferable or 
optimal linear scales for SSF search are L=1.5 Mpc for HeH + (I=0) at z=20 
and L=0.75 Mpc for LiH(I=6) at z=30, correspondent angular sizes according 
to (7) are 10.7' and 7.5'. 

Recent WMAP polarization measurements have shown strong probability of 
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Table 5 

Spatial scales L of peaks in CMBR power spectrum for different z 
z Li, Mpc L 2 ,Mpc L 3 ,Mpc L 4 ,Mpc 



1000 


0.22 


0.11 


0.073 


0.043 


100 


1.99 


0.99 


0.66 


0.39 


50 


3.72 


1.86 


1.24 


0.74 


30 


5.79 


2.91 


1.93 


1.16 


20 


8.08 


4.04 


2.69 


1.62 



an early reionization scenario with z r = 20(+10/ — 9) (Kogut et al, 2003, 
Page et al, 2006) due to Pop III stars (Wyithe, Loeb,2003) in Galaxy halos. 
SSF detection at z=20-30 may provide us with important information about 
ionization history of the Universe. 

The question is which maximum proto-object scales are reasonable, could exist 
at z=20-30 or are connected with some known scales in the Universe. Let us 
compare them with spatial scales of primary CMBA. The angular scale of the 
first peak in primary CMBR power spectrum 9i ~ 50', the following CMBA 
peaks are ~ l/2#i,l/3#i, and l/5#i (see Fig. 2 down). Then possible spatial 
scales L of 1-4 peaks in CMBR power spectrum for different z are given in the 
Table 5 with (7) and (6). 

The above consideration shows that at AT/T=2-10 ~ 5 level and with 
Av/is=2% we can search for SSF which are the result of large scale proto- 
object motion with peculiar velocity relative to CMBR at z=20-30. These 
possible proto-objects are still less than spatial scale of the fourth peak in 
CMBR power spectrum (Table 5) which might be related to a proto-galaxy 
scale. With less Av/is, higher sensitivity and multipole resolution we can ap- 
proach to smaller spatial scales of proto-objects up to proto-star scales. 

SSF firstly proposed as a new phenomena in (Dubrovich, 1977) is one of 
most promising secondary CMBA. SSF is non-isotropic spectral distortions 
of CMBR which are to be present at z=300-10. Homogeneous and isotropic 
spectral distortions of CMBR may be present at z=6000-1000 due to pri- 
mordial Helium and Hydrogen recombination (Dubrovich, 1975, Dubrovich, 
Stolyarov, 1995, Dubrovich, Shahvorostova, 2004, Dubrovich, Grachev, 2004). 
Non-isotropic distortions of CMBR similar to SSF may be present in this 
epoch at lower AT/T level as well(Dubrovich, 1977/2). Calculation of this 
effect in more details at z=1000-1100 is done in (Rubino-Martin, Mernandez- 
Monteagudo and Sunyaev, 2005). 
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3 Simulation of SSF 



Let us remind that CMBR temperature fluctuations are expanded in spherical 
harmonics (Peebles, 1993, Peacock, 1999): 

oo I 

AT(#,0)/To = £ E o. lm Y lm {9,4>), (10) 

1=0 m=-l 

The angular power spectrum Ci of AT/T can be related to the autocorrelation 
function C{6) as: 

C{6) = 1/4^(2/ + l)C,P,(cos0), (11) 
i 



where P; are Legendre polynomials, 1-multipole moment and C\ —< \ai m \ 2 >. 
Theoretical predictions of AT are usually quoted as 

(AT,) 2 = 1(1 + l)Q/2n, (12) 

AT = + 1)Q/2tt, (13) 



As we consider small-angle patches of the sky CMBR fluctuations AT = T— < 
T > may be generated by evaluating the simple Fourier series (Bond et al, 
1987): 

n u -in v -i 2-jr 
AT(9 x ,9 y )/T = J2 E D(n u ,n v )exp[t—(n u 6 x + n v e y )] (14) 

n u =0 n v =0 



where L denotes the size of the simulated region in radians; (6 x ,0 y ) are 
Cartesian coordinates on the sky (spatial domain); (n u ,n v ) are coordinates 
of Fourier components D in spatial frequency domain, < \D(n u , n v )\ 2 >= Ci, 
l = 2n/L-^/nl + nl 

For the search of SSF CMBR maps obtained in several spectral channels are 
to be studied and then processed by the special frequency- differential or the 
difference method (Dubrovich, Bajkova,2004). The Difference method can be 
considered as an alternative method to the correlation function analysis. In this 
case we are interested in the first derivative of spatial distortions with respect 
to the frequency. The method used is based on the analysis of a difference of 
two CMBR temperature maps observed at different frequencies and reduced 
to one beam shape. Mainly, such a difference map contains information only 
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on the secondary fluctuations because the primary CMBR fluctuations present 
in both maps will be eliminated due to their black-body spectrum nature. 

The procedure of SSF simulation is the following. Firstly CMB map is gener- 
ated with a standard CMB power spectrum. Then amplitudes of FFT spec- 
trum of CMB map are kept while phases are modified as random values uni- 
formly distributed in the interval 0-27T. A backward FFT gives us an approx- 
imation of CMB map which contains primary and secondary fluctuations in 
conditions of unknown theoretical power spectrum of SSF. The actual theo- 
retical power spectrum of SSF is to be a subject of nearest energetic efforts, 
it is strongly needed for following SSF investigation. 

Let the difference of the CMBR observation frequencies be equal to Aui. Let 
the limiting angular fluctuation size 6\ corresponds to this frequency difference 
in accordance with equation (5). Obviously, the fluctuations of size larger than 
Q\ will be seen at both maps. After subtraction of one map from the other 
the fluctuations of size larger than Q\ will be mutually suppressed and the 
remainder map will contain only the fluctuations smaller than 6\. Evidently, 
this fact can be seen also from angular power spectrum of this map which is 
close to zero at frequencies with multipoles I < l-y and different from zero at 
I > h, where l\ m 1/Q\ is the break of the spectrum. Now let the difference of 
two frequencies be A 1/2, where Az/2 > Av\. The angular size 82 corresponds to 
this frequency interval. Obviously, 62 > Q\ in accordance with equation (5). In 
this case we will see at the difference map the fluctuations with size smaller 
than 6 2 and the corresponding angular power spectrum will have the break at 
l 2 « l/6 2 < h- 

Our CMBA model also includes convolution of expected CMBA maps with 
a simulated radio telescope beam and adding a pixel noise. Differential maps 
are preliminary processed to decrease the noise in given power spectrum re- 
gion. Let convolved CMBR maps have spatial scales \ min -r- (\ max -mAl m ) , then 
the spatial scale of differential maps is (l^^-nA^)^^^, where m=0,l,2..., 
Al m /1 corresponds to AA m /A(2% in our case). Let us limit l max in the sim- 
ulation by 1=1056 which is near to the lowest reasonable limit of \ max in the 
ground based MM-wave SSF experiment. We will use in simulation a stan- 
dard ACDM CMBR model (Spergel et al., 2003, Spergel et al, 2006) with 
parameters: Q b x h 2 =0.2(H /l00 km/sMpc, H = 65 km/sMpc), fi A =0.65, 
Qm = 0.3, n=l, h=0.65. Initial single frequency CMBR maps (Fig. 2 up) 
15°xl5° (128x128 pixels) have 1=3-1056(1), 1=3-1032 (2), 1=3-1008 (3), 1=3- 
984(4), 1=3-960(5), 1=3-936(6), 1=3-912(7), all ^(quadratic sum)=25902. Dif- 
ference maps (Fig.2 down) have 1=1032-1056 (1-2), (7=14.4, 1=1008-1056(1-3), 
(7=21.3, 1=984-1056(1-4), <r=32.4, 1=960-1056(1-5), a=42.9, 1=936-1056 (1-6), 
cr=60.8, 1=912-1056, <r=87.6 with the angular power spectrum corresponding 
to ACDM model for initial (solid line) and differential (dashed line) CMBR 
maps (Fig.2 middle). 
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Fig. 2. Initial single frequency CMBR maps 15°xl5° (128x128 pixels) (up), dif- 
ference maps (middle), power spectrum of AT corresponding to ACDM model for 
initial (solid line) and differential (dashed line) CMBR maps (down) 
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It was shown in (Dubrovich, Bajkova,2004) that the Difference method may 
be used when S/N ratio for each map is 1.44. In case of higher level of noise 
a more complicated preliminary signal processing directed to decreasing the 
noise in the given CMBR spectrum region would be required. Fig. 3 demon- 
strates predicted difference SSF map without an instrumental noise(up) and 
with white instrumental noise (down) added in single frequency maps with 
S/N=1.17. Power equalization method of Wiener filtration (Gorski, 1997) was 
applied in this case as signal preliminary processing to decrease the noise in 
given power spectrum region. In spite of a correlation coefficient of two dif- 
ference maps of Fig. 3 is less than 0.15 the power spectra of the difference 
signal after Wiener filtration are similar (Fig.3 right) that is the most impor- 
tant to recognize SSF in a real SSF search experiment. But this is true for 
white instrumental noise and we concentrate our efforts in section 4 to avoid 
or minimize abnormal instrumental pixel noise in the SSF search experiment. 

It was assumed in above simulation that the observed CMBR maps consist 
only of the primary and secondary CMBR fluctuations with some additive 
white instrumental noise. However, in practice the CMBR detected maps in 
addition contain a number of astrophysical foregrounds such as extragalactic 
unresolved sources, free-free and synchrotron radiation, SZ effect, Galactic 
dust and so on. Astrophysical foregrounds are characterized by considerably 
low value of spectral index as compared with the analyzed secondary CMBR 
fluctuations and are successfully eliminated in difference maps. 



4 A concept and strategy of SSF Experiment 

Up to now practically all CMBR ground based and space experiments study 
primary CMBA which has black body spectrum. Primary CMBA was 30 years 
under investigation at RATAN-600 radio telescope since 1977 (Parijskij, Ko- 
rolkov, 1986, Parijskij et al, 2005). 

First attempts of SSF observations were done in 1992 at IRAM MM radio tele- 
scope with the upper AT/T limit 2-10~ 3 at 1.3 mm (De Bernardis, Dubrovich 
et al, 1993) and in 2001 at RATAN-600 (the near pole survey) with the upper 
AT/T limit 10 -3 at 6 cm (Gosachinskii, Dubrovich et al.,2001). Surplus fre- 
quency resolution of the existing spectrum-analyzers (30 MHz and 20 KHz) 
increased AT/T limit in both cases. The Far Infrared Absolute Spectrome- 
ter(FIRAS) instrument has shown that CMB has ideal black body spectrum 
from 60 GHz to 600 GHz and CMB temperature is 2.725+-0.001 K(Fixsen, 
Mather, 2002). Then Fixsen et all reported about an effort to measure dis- 
tortions of CMB black body spectrum due to of various processes in Early 
Universe. This attempt was done in ARCADE balloon-born experiment where 
CMB temperatures were measured at 10 GHz and 30 GHz with accuracy +- 
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Fig. 3. Predicted difference SSF map without noise(left up) and with white noise 
(left down), S/N=1.17, the correspondent power spectrum without noise - solid line 
and with white instrumental noise - dashed line (right) 




0.010 K and +-0.032 K (Fixsen et all, 2004). In this paper we suggest a concept 
and strategy of the SSF experiment with AT/T limit 2 • 10~ 5 . High absolute 
measurement accuracy is not necessary in this experiment as we need to mea- 
sure differential temperatures in nearby receiver channels without measuring 
absolute CMB temperatures. 
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Fig. 4. Multibeam Solar Radio Telescope (MSRT) of Tuorla Observatory, Finland 




Simulation of possible SSF experiment is done for MSRT (Khaikin et al, 
2002,2003), Tuorla Observatory, Finland (Fig. 4) equipped with a 7x4 beam 
cryo-microblometer array at 3 mm with RF multiplexer providing up to 7 
spectral channels in each beam. The scheme of the experiment is given in 
Fig. 5. Firstly idea of such an experiment is suggested in (Dubrovich, Bajkova, 
Khaikin, 2004, Khaikin, Dubrovich, 2004, Khaikin, Luukanen, Dubrovich, 
2005). In this paper we develop the concept and present simulation results 
of a possible experiment. 

RF filters 

The challenge is to split RF band in front of bolometers. If superconducting 
Nb microstrip resonant A/4 stub filters are used in RF multiplexer (RFM) 
the in-band transmission can approach unity due to negligible loss of the 
superconducter. To simplify a multiplexer design the number of filters and 
spectral channels may be reduced to 5 or even 3. A wideband microstrip 
superconducting frequency multiplexer is applied in front of TES bolometers 
in SAMBA(Superconducting Antenna-coupled Multifrequency Bolometric Ar- 
ray), where low loss and excellent out-of-band rejection is demonstrated at 3 
frequencies of MM band (Goldin et al, 2003). Antenna-coupled TES bolometer 
with an integrated microstrip BP filter at 217 GHz is applied for a bolometer 
array in a ground based CMB polarization experiment (Myers et al, 2005). 
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Another way to split radio frequency band in a bolometric receiver is to build 
a quasi-optical filter based on metal mesh filters(MMF) or frequency-selective 
surfaces (FSS). A review of quasi-optical MMFs is given in (Ade et al, 2006). A 
few possible designs of cryogenic MMFs with bandwidth up to 6% at 1-3 mm 
are given in (Page, 1994). 3 cryogenic filter wheels with several bandpass filters, 
are under development now in Columbia University for Millimeter Bolometric 
Array Camera (MBAC) to provide 5-10% frequency resolution in one octave 
MM band from 150 GHz to 300 GHz at Atacama Cosmology Telescope(ACT). 

FSS may be also used to separate in space radio telescope foci in which several 
feeds or arrays at different frequencies may be placed (Kasyanov, Khaikin, 
2004). But it is difficult to achieve frequency resolution with MMF or FSS 
better than 5% and superconducting N b microstrip resonant filters are more 
effective in our case. 

Array feeds 

We consider several variants of an array feed. As a narrow flare corrugated 
horn is difficult to build in 3 mm band a trimode conical horn may be used 
to provide a small flare angle with an acceptable sidelobe level (less than -20 
dB in E/H plane) and low frequency cut-off near 87 GHz. A smooth- walled 
spline-profile horn(Granet et al, 2004) is an alternative to both corrugated 
and trimode horns and it allows us to achieve an identical radiation pattern in 
E/H plane in the band, low sidelobe level(less than -25 dB), more tight array 
packaging and lower manufacture cost. We are in final stage now to build and 
test such a horn for array applications al lower frequencies (30 GHz-38 GHz) 
in cooperation with C.Granet(CSIRO) (Khaikin et al,2007). 

Multipole resolution 

A 2 m precise carbon fiber MSRT dish (Khaikin et al,2002) allows us to work 
in 100 GHz and 150 GHz bands but the second one is much more difficult 
due to higher atmospheric attenuation and antenna scattering background. 
Angular resolution of MSRT is 7'-& at 88 GHz-100 GHz that is enough for 
a large scale SSF experiment and corresponds to l ma:r =7r/#=1542-1800 (for 
full sky coverage) the same as at Plank (7', 1=1542) and twice better than 
at WMAP(l ma:r =750). Nevertheless for a ground based or balloon experiment 
Al < I < l m ax — AZ, where Al = In jM (size of the symmetrical map) and 
actual spatial resolution of MSRT is not higher than 1=1200-1500 for partial 
sky coverage. In Boomerang balloon experiment actual resolution was less 
than 1=1025 at 150 GHz(de Benardis et al.,2002). 1000> I < 1500 seems to 
be optimal for SSF search at current stage. At higher 1 with the low frequency 
resolution S/N ratio falls. Besides, in ground based experiments with higher 
angular resolution such as Cosmic Microwave Imager (CBI) (Readhead et 
al,2004) with 1 rnax = 3500 point sources were a serious problem at 26-36 GHz 
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as the power spectrum of primary CMBA falls rapidly with 1 while the point 
source contribution grows as l 2 in the band. Point source contamination is still 
dangerous at 3 mm at higher 1 if secondary CMBA has a falling spectrum. 
Some secondary CMBA may have inverted C/ spectra when the magnitude of 
Ci grows at higher 1 (Doroshkevich, Dubrovich, 2001). However SSF search at 
1>1500 needs both higher sensitivity and higher frequency resolution which 
contradicts each other. 

Chopping mode 

A multifeed or receiver array is needed in a SSF experiment to get a CMBR 
map in an unmovable radio telescope mode (preferable) or to accelerate map- 
ping in a scanning mode (less preferable because of variable instrumental ef- 
fects). We plan to use a chopping mode to reduce an abnormal instrumental 
noise and atmospheric effects. A chopper with beam deviation 9 c h will distort 
power spectrum of SSF due to filtration of scales 9 > 9 c h and changing the 
amplitudes at scales 9 < 9 ch but this is not so dangerous in the task of SSF 
search. The main role of the chopper is to suppress uncorrelated residual 1/f 
noise of different bolometers which can not be eliminated in difference maps. 

Calibration 

Identical and stable calibration of all receiver channels is very important in this 
experiment. It may be produced quasi-optically by a lens, forming a lateral 
plane wave for a whole array. A high stable noise generator with a direct 
coupler or a matched black body load is placed in the lens focus. A calibration 
source must be put in a reverse two level thermostat realizing 0.1 K accuracy 
at the first and not less than 0.01 K accuracy at the second level. 

CMBR map 

According to our concept MSRT central beam is fixed near the North Celestial 
Pole (NCP) at 5=86° to increase cosmic source passing time via the beam up 
to 348s -402s and to observe the 1° x 8° strip preferably in low culmination 
from a=8 /l 00 m to 16 A 00 m at high Galactic latitudes. Actual spatial resolution 
of a combined CMBR map in the scheme of Fig. 5 is 6'-7'(HPBW). The long- 
length focus Cassegrain scheme with m=(e+l)/(e-l)=10-ll (eccentricity e>l) 
gives us negligible aberrations of off-axis beams but requires more complicated 
array feeds. A short-length focus Cassegrain radio telescope scheme gives us 
unacceptable level of off-axis aberrations(Khaikin, Luukanen, 2003). 

Antenna characteristics 

We simulated some antenna characteristics for the scheme of the experiment 
given at Fig. 5. Antenna parameters are D=2 m, d=0.3 m, F/D=0.435, interfo- 
cal distance f=150cm, feed flare angle is 12 deg. Simulated characteristics are 
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Fig. 5. Suggested scheme of the experiment (up), multichannel bolometer receiver 
(MBR) scheme with RF multiplexer (down) 
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given in Table 6. Multibeam Focal Plane Array (FPA) configuration presented 
in Fig.5 allows us to avoid undersampling in a survey mode. For simulation of 
the multibeam radiation pattern, calculation of some other antenna character- 
istics (illumination, aperture, spillover, sidelobe efficiencies, beam overlapping 
level etc.) FOcal Plane Array Simulation (FOPAS) program package has been 
used (Khaikin, 2005). FOPAS takes into account geometry of a single or dual 
reflector radio telescope, expected aperture distribution or the beam pattern 
of an ideal broadband feed, the feed removal from the focus and off-axis aber- 
rations. FOPAS allows us to optimize FPA geometry with the beam spacing 
and the overlapping level which are optimal for the given astrophysical task. 
The task of multibeam FPA optimization is also considered in (Majorova, 
Khaikin, 2005). Multibeam radiation pattern simulated with FOPAS is given 
in Fig. 6. 



TES bolometers 



We consider Transition Edge Sensors (TES) as cryo-bolometers for this exper- 
iment (Khaikin, Luukanen, Dubrovich, 2005). The lateral dimension of TES is 
less than 10 microns (microbolometer) and standard lithographic techniques 
may be used to build arrays of 100-1000 sensitive elements. 

Photon and phonon Noise Equivalent Power (NEP) and Noise Equivalent 
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Table 6 

Simulated antenna characteristics 



Number 
of beams 


Feed 

taper, 

dB 


HPBW at 
88-100 GHz, 
arcmin 


FSL, 
Db 


Beam 

overlapping 

level, 

dB 


28 


-10/-13 


7-6 


21/24 


-0.45 


Spillover 

efficiency, 

dB 


Blockage 
efficiency 
dB 


Illumina- 
tion effi- 
ciency 


Ruze effi- 
ciency (94 
GHz) 


Antenna ef- 
ficiency (94 
GHz) 


-0.45 


-0.23 


0.91 


0.75 


0.51 



Fig. 6. Simulated multibeam radiation pattern of MSRT at 100 GHz, the feed taper 
is -13 dB 




Temperature Difference (NETD) of TES for the ground based SSF experi- 
ment are considered in (Khaikin, Luukanen, Dubrovich, 2005). An expected 
phonon NEP of a cryo-microbolometer with the physical temperature of 0.3 
K is <1T0~ 17 W/Hz 1//2 . We do need the physical temperature near 0.3 K 
(Fig. 5) as NEP of TES dramatically grows at higher physical temperatures 
(Luukanen, Pekola,2003). Photon NEP from radio telescope antenna and at- 
mosphere is not less than 2-10 -17 in good atmospheric conditions. This will 
give us o"=1.0-2.0 mK s 1 / 2 sensitivity per beam per channel (2 GHz band- 
width) in the best atmospheric conditions with T atm <20 K (the lowest atmo- 
spheric opacity). In 8x8 element planar bolometer array at 90 GHz for GBT 
operating temperature is 450 mK, which limits NEP of TES bolometers to 
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2xl0~ 17 W/Hz 1 / 2 (Benford et.al,2004). The physical temperature T p of an 
antenna coupled with a superconducting bolometer must be chosen to provide 
a photon NEP<10~ 17 W/Hz 1//2 . For rather high radiative antenna efficiency 
77=0.9 and T p =300 K an additional photon noise from an antenna is 10~ 16 
Wt/Hz 1//2 and cooling of antenna (horn)up to at least 10 K is necessary to 
provide its own photon NEP at less than 10~ 17 Wt/Hz 1/2 level (Khaikin, 
Luukanen, Dubrovich,2004)(Fig.5). For deep cryogenic cooling Sorption He 3 
coolers at 0.3 K of Applied Physics Institute of RAS in Nizhnij Novgorodmay 
be used(Vdovin, 2005). 

Some ground based CMBR projects (Lawrence et al,2004) chose an alternative 
way and are being built cryogenic MMIC MM-wave cameras with expected 
sensitivity of 500 jiKs 1 / 2 in continuum in 3 mm band. We have considered such 
alternative as well. InP MMIC LNA at 15 K with a broadband subharmonic 
mixer and IF microstrip 10 channel multiplexer may give us 4-5 mKs 1 / 2 per 
beam per channel(l% bandwidth) at 85 GHz -105 GHz (Khaikin et al, 2004). 

Sensitivity in the experiment 



The actual reachable integration time of a bolometer depends on the Knee 
frequency f& n in its noise power spectrum. In a bolometer case f& n may be 
reduced from tens of Hz to tens of mHz by biasing the bolometer with an 
AC current and synchronous demodulation in readout electronics (Crill et al, 
2002). ffc =10 mHz was reached by AC biasing of bolometers in a 144 pixel 
cryo-bolometer array BOLOCAM at 140-270 GHz (Mauskopf et al., 2000). 
With AC biasing of the TES bolometers we expect ffc n =20 mHz. In the effective 
chopping mode with chopper frequency f c /i=5 Hz ffc„=5 mHz may be reached 
then 3-4 minutes integration time becomes available. Using a 100 ms sampling 
rate for f c /=5 Hz and 170 s synchronous integration time we can reach a of 
thermal noise not worse than 0.1-0.2 mK per beam per pixel. Four l°x8° maps 
(128x7 pixels) are to be recorded during about 8 hours in clear cold winter 
nights (after Sunset and before Sunrise) and then combined into one map 
(64x14 pixels) with minimum angular resolution 6' x 6' and a noise factor of 
1/2 1 / 2 in comparison with initial maps. 20-25 identical maps must be summed 
for 20-25 clear cold nights with low atmospheric opacity in order to reach 
30/iK a of thermal noise in the averaged single frequency maps. We expect 
contribution of non thermal instrumental noise, atmosphere and foregrounds 
in differential maps up to 30 /i K as well. So we may expect total a =50-55 
fiK in differential maps. 

Atmospheric limitations 



Big efforts have to be made to reduce the contribution of instrumental and 
foreground effects at AT/T=10~ 5 level that corresponds to 10~ 7 of ambient 
temperature. Some instrumental effects such as 1/f noise, ground spillover, ra- 
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diometric offset may be comparable with an atmospheric noise but their con- 
tribution may be significantly reduced in the ground based experiment. Atmo- 
spheric fluctuations are to be effectively suppressed with a chopper while high 
atmospheric attenuation is the most dangerous for a ground based MM-wave 
SSF experiment with cryo-bolometers. Atmosphere gives the main limitation 
in sensitivity in MM wave ground based SSF experiment with cryo-bolometers 
as it will contribute up to 80% in a system temperature and in some cases 
may saturate (overheat) a cryo-bolometer. The bolometer must be able to ab- 
sorb up to 1 pW in order to avoid saturation with photon NEP > 2 x 1CT 17 
Wt/Hz 1 / 2 . This requires a rather high dynamic range D which is the ratio of 
the maximum power detected to the minimal detectable power. For expected 
sensitivity 1 mK s 1 / 2 D must be not worse than 45 dB. In fact only days with 
lowest atmospheric opacity (lowest water vapour content) may be used for 
such an experiment at sea level. It was shown by a "tipping curve" method 
that typical atmospheric opacity at MSRT site in zenith r z is 0.2 dB at 36.5 
GHz for dry clear Autumn nights (Khaikin et al, 2003) that gives predicted 
r 2 at 94 GHz 0.3 dB - 0.5 dB in clear cold winter nights with lowest inte- 
grated water wapour content (< 2g/m 3 ). An effective beam wobbling made 
it possible to reduce residual atmospheric noise AT atm and spillover AT sp m 
up to 10 fiK at 40 GHz in SK CMBR experiment in Canada (Wollack et al, 
1994), where the radiometric offset was less than 7 /iK per day for all obser- 
vations. A special ground screen (height up to 3 m) is to be installed near 
MSRT (Fig. 5) to avoid the contribution of the hot atmosphere and ground 
via sidelobes, spillover and scattering background and provide Tao= 33 K at 
high elevation angles including own antenna temperature, atmosphere (20 K) 
and CMBR(Khaikin, Luukanen, Dubrovich,2005). The residual non-thermal 
instrumental noise must be significantly reduced up to 30 /iK in differential 
maps while thermal noise grows by a factor 2 1 / 2 in comparison with single fre- 
quency maps. Non-thermal noise can not be reduced lower in difference maps 
due to residual abnormal (1/f) instrumental noise, some difference in radia- 
tion patterns, in amplitude-frequency characteristics of the channels and some 
difference in atmospheric effects within 12% bandwidth. 

In spite of atmospheric limitations a ground based SSF MM wave experiment 
has advantages in the cost, flexibility and possible duration. In space condi- 
tions TES bolometers may give one - two orders better sensitivity. In (Ali et 
al, 2003) NETD=25 /i K/s 1/2 per pixel is expected in space for CMBR at 90 
GHz for 100 TES detectors with the physical temperature less than 0.1 K 
and NEP=2 x 10~ 18 Wt/Hz 1 / 2 . Yet, most of new generation MM wave CMBR 
projects chose multi-element bolometer arrays to reach highest sensitivity in 
the ground based experiment. NET C mbr = 260 — 280/ifTs 1 / 2 per detector 
must be reached in 1000 element array of TES bolometers at 10 m South 
Pole Telescope(SPT) in 20% bandwidth at 95 GHz and 150 GHz(Ruhl et al, 
2004). Millimeter Bolometric Array Camera(MBAC) with 3 1000-element TES 
bolometer arrays at 150 GHz/220 GHz/270 GHz is under development now for 
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Table 7 

Expected characteristics of the experiment 



Alpha, 
hours 


Delta, 
degrees 


Multipol 
moment 1 


Frequency 

band, 

GHz 


Frequency 
resolution, 

% 


8-16 


86 


1200-1500 


88-100 


2 


Sensitivity 
per beam, 
mKs 1 ^ 


Integration 

time, 

s 


Integral sen- 
sitivity of the 
experiment, 
micro K 


dT/T 


Duration of 
the experi- 
ment 


1-2 


170 


50-55 


2xl0~ 5 


Winter sea- 
son 



6 m Atacama Cosmology Telescope(ACT) Project (Fowler, 2004). Expected 
MBAC sensitivity to CMBR is 300 /iK/SOO^K/TOO^Ks 1 / 2 per detector 

Astrophysical Foregrounds 

In a ground based experiment at sea level atmospheric effects contribute much 
more than astrophysical foregrounds nevertheless we should take them into 
account to better prepare the experiment. Galactic and extragalactic dust 
emission will significantly dominate at 3 mm over other foreground effects such 
as discrete and diffuse sources, galactic synchrotron and free-free emission, hot 
clusters of galaxies. Dust fluctuations are seen on all angular scales from 2' 
to tens of degrees on IRAS data (Low et al.,1984). Antenna temperature is 
T a ~ v l3d , where Pa is the spectral index which depends upon the emissivity of 
the dust and lies in the range of 1.5-2.2 according to COBE, IRAS and WMAP 
data. After 60 GHz Galactic and extragalactic dust dominate and we should 
carefully chose CMBR map coordinates in dust free areas of the Universe at 
high Galactic latitudes. The full sky thermal dust map based on data from 
IRAS and COBE using MEM procedure shows a wide range of temperatures 
from 400 fxK to 15 fiK in W band of WMAP (94 GHz) (Bennet at el, 2004). 
Contribution of dust emission is to be significantly eliminated in differential 
maps as this component must be well enough correlated in single frequency 
maps due to its black-body spectrum. 



5 Final remarks 

As it was shown in section 2 the expected signal (AT/T of SSF) depends 
on many Cosmological parameters and conditions in Early Universe such as 
fl m ,flA, Qb, n Q H,®,Vp and others while the best S/N ratio at given z also de- 
pends on Au/u. It is extremely difficult to predict the most probable combina- 
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tion of these parameters at different z in Early Universe and SSF experiments 
are urgently needed to set SSF an uppper limit at AT/T=10 -4 — 10~ 5 level 
and specify theoretical predictions of SSF. 

We want to emphasize in final remarks the significance of SSF detection which 
leads to the discovery of the Cosmological molecules. Cosmological molecules 
can provide an unique test to non-standard nuclear synthesis at high z. SSF 
can only help us in investigating Dark Ages epoch in the Early Universe (300 > 
z > 10). SSF detection at z=20-30 may give an important information about 
ionization history of the Universe. Finally, they will allow us to explore physical 
properties of Dark Matter and Dark Energy. 



6 Conclusion 

Observational effects caused by interaction of primordial molecules with 
CMBR seem to be most promising in investigating of Dark Ages epoch of 
the Early Universe and ionization history of the Universe. Strategy and con- 
cept of the experiment are given. Simulation of the experiment in an atmo- 
spheric transparency window at 3 mm with MSRT equipped with a 7x4 beam 
cryo-microbolometer array gives us the expected sensitivity of 50-55 fiK in 
the differential CMBR maps near NCP with angular resolution 6'- 7' (1=1200- 
1500) and frequency resolution Av/u—2%. The expected integral sensitivity 
in the experiment is AT/T=2-10 ~ 5 which is near to optimistic estimates of 
SSF caused by molecules HeH + at z=20-30. Simulation shows that SSF may 
be recognized in the angular power spectrum when S/N in single frequency 
CMBR maps is as small as 1.17 for a white noise. Such an experiment will 
give us a possibility to set upper limits in SSF observations in MM band and 
prepare future SSF experiments. 
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